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Abstract 


In this work, the tv\o dimensional a\is>mmetric. incompressible steads flow and 
forced com ectis e heat transfer of a non Newtonian pow er-law fluid across a bank of circular 
cylinders has been investigated numerically. The inter-cylinder interactions are modelled 
using the Happel's free surface model. The equations governing fluid flow and convective 
heat transfer are solved using the finite difTerence method. Simulations, to obtain the values of 
the local and average Nusselt number, were carried out over a wide range of phvsical and 
kinematic parameters as :0,4()< jc <0.60. 0..5<n<l. l<Re<500 and l^Prs .SOO for both 
boundary conditions, i.e., the constant wall temperature and the constant heat flux at the 
surface of the cylinder. The angular variation of the local Nusselt number on the surface of the 
cylinder as well as isotherms are shown under vary ing conditions of porosity. iK>vver-lavv 
index. Pnmdtl and Reynolds number to gain irseful insights in to the phy sics of the flow and 
heat transfer. The results obtained agree well with the results for limiting conditions available 
in the literature. Limited comparisons with experimental data are also presented. 
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Chapter 1 


INTRODUCTION AND SCOPE OF THE WORK 

1.1 Introduction 

The steady, axisymmetric and incompressible flow of non-Newtonian fluids 
nonnai to an array of cylinders and heat transfer from a collection of cylinders 
immersed in moving fluid streams represents an idealisation of many industrially 
important processes encountered in chemical, mineral and process engineering 
applications. Typical examples include flow in packed and fluidised beds 
employed for drying of fibrous materials (coconut shell, rice husk, etc.), in the 
shell of tubular heat exchangers, in aerosol and cigarette filters and during the 
transport of sap in plants and trees (Tripathi and Chhabra, 1996). Other examples 
are to be found in polymer process engineering applications such as the 
manufacture of fiber -reinforced composites entails the flow of a resin through a 
bed of fibers (Willians, Morris, &Ennis, 1974; Gupta, 2000). There is no question 
that the variable of central interest in all these applications is the frictional 
pressure loss and the rate of heat transfer for a given flow rate and the thermo - 
physical properties of the fluids. Over the years, several anal 3 d:ical treatments 
have been performed for the creeping flow of a incompressible Newtonian fluids 
normal to a bank of cylinders arranged in a variety of geometric configuration, 
eg., triangular, square, etc. For finite values of the Reynolds number of flow, the 
governing equations become non-linear and only limited numerical results are 
available for pressure drop and heat transfer in the intermediate Reynolds number 
region (Reynolds number up to 100). Limited comparisons suggest that the 
correspondence between predicted and experimental data is satisfactory for the 
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flow of Newtonian fluids (Skartsis, 1992). Most of the theoretical developments 
in this area have been summarised by Drummond and Tahir (1984) whereas the 
pertinent experimental studies have been reviewed by Skartsis et al (1992).The 
corresponding studies on heat transfer have been dealt with by Zakuskas (1987) 
and by Krieth (2000). 

The heat transfer from individual circular cylinders oriented normal to the 
flow has also been the subject of analytical investigations, but these studies have 
not been as intensive as those involving the flow phenomenon especially for non- 
Newtonian fluids. The recent work of Mandhani et a/.(2002) on the heat transfer 
from an array of circular cylinders in cross flow was the recent numerically 
simulated in the range of Reynolds number (1-500) and prandtl numbers (0.7 and 
7.70) for Newtonian fluids with constant themio-physical properties, is first 
extensive study. And also recently Bailer et al. (1997) investigated numerically 
heat transfer and air flow around a single heated cylinder for Reynolds number of 
200 and 2000. Similarly, Lange et al. (1998) have studied numerically the 
momentum and heat transfer characteristics of air flow around a single cylinder. 
Most of these studies have been reviewed recently by Ahmad (1996). From a 
theoretical standpoint, Sangani and Acrivos (1982) calculated the average 
temperature difference between the bulk fluid and the uniformly heated cylinders 
under conditions of small Reynolds number (Re=0) and Peclet numbers. Hsu 
(1964) analysed the phenomenon of heat transfer to liquid metals (high Prandtl 
number) in cross flow through rod bundles by invoking the inviscid flow 
approximation, as under these conditions, the momentum boundary layer would 
be much thinner than the thermal boundary layer. 

In contrast to this, little is known about the flow of non-Newtonian 
fluids over tube bundles such as that encountered in the production of fiber 
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reinforced resins using an autoclave process, in heat exchangers employed to 
heat/cool non-Newtonian fluids in food and biotechnological applications and in 
enhanced oil recovery operations. Furthermore, this flow configuration has been 
used to capture some aspects of non-Newtonian fluid flow (such as extensional 
effects) in porous media. Despite such overwhelming theoretical and pragmatic 
significance, little research has been directed at the elucidation of the role of non- 
Newtonian characteristics in this flow configuration. Metzner et al. (1960) have 
investigated the heat transfer to non-Newtonian fluid systems passing though 
smooth round tubes under laminar conditions. Metzner and Reed (1955) have 
derived a generalized Reynolds number to correlate the friction factor for power 
law non-Newtonian fluids in cylindrical conduits. Metzner et a/.(1957) presented 
theoretical analyses combined with an experimental study of the variables 
determining the heat transfer coefficient for laminar flow of non-Newtonian fluids 
in pipes. Shah et al. (1962) obtained analytical solutions and experimental data for 
free convection heat transfer from a single cylinder to a power law non- 
Newtonian fluid. The purpose of the work was to determine the feasibility of 
application of rheological equations of state derived from viscometer data to 
predict the fluid behaviour in laminar boundary- layer flows. Adams and Bell 
(1968) investigated the fluid friction and heat transfer for flow of CMC solutions 
across banks of tubes. 

From a theoretical stand point, a mathematical description of inter-cylinder 
interaction is also required in addition to the usual mass, momentum and energy 
conservation equations. In the literature, there have been two approaches to model 
the strong particle interactions. The first approach attempts to solve the field 
equations for a preconceived geometrical configuration such as triangular or 
square array in-line or staggered configuration such as that employed by Sangani 
and Acrivos (1982). Notwithstanding the significance of detailed kinematics, it is 
generally believed that the macroscopic quantities such as drag on cylinders/ or 
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pressure drop in fibrous filters and the overall rate of heat transfer are relatively 
insensitive to the actual arrangement of cylinders and are primarily determined by 
the value of the porosity of the system and the kinematics conditions (namely, the 
Reynolds number Prandtl number, porosity, rheology of the flow). Thus the cell 
model approach which converts a difficult many-body problem into a 
conceptually much simpler one-body equivalent thereby making it more attractive 
both analytically as well as numerically affords a reasonable analysis of such 
extremely complicated flows. This approach will be used in this work to study 
heat transfer characteristics between a tube bundle and a fluid stream for non- 
Newtonain fluids. 

1.2 Scope of the work 

This thesis sets out to solve the mass, momentum and energy equations for the 
incompressible steady flow and energy equation of non-Newtonian fluids in a 
collection of cylinders via the so called free surface model. In particular the 
equations of continuity and momentum together with appropriate boundary 
conditions have been solved numerically to obtain the detailed kinematics of the 
flows. In this study, the equation of energy for the Reynolds number (1-500) and 
Prandtl numbers (1-500) for non-Newtonian fluids with constant thermo physical 
properties, have been solved numerically. Extensive results on detailed 
temperature profiles as well as on the local Nusselt number are obtained over 
wide ranges of physical (Prandtl number, porosity of arrays) and kinematic 
(Reynolds number) conditions. Furthermore, the role of the usual two types of 
thermal boundary conditions on the cylinder, constant temperature and constant 
heat flux, has been elucidated in determining the values of the Nusselt number. 
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Chapter 2 


PROBLEM STATEMENT AND GOVERNING 

EQUATIONS 


In this chapter, the flow problem investigated herein is outlined briefly, 
followed by the model development. This then leads on to the development of 
the field equations, namely, the momentum, energy and continuity equations 
together with the specification of physically possible boundary conditions 
employed in this study. 

2.1 Problem Statement 

The steady incompressible flow of a power law fluid past an array of 
long circulsff cylinders, which are oriented normal to the direction of flow, as 
shown schematically in Fig (2.1a) is considered. The rheological equation of 
state for a fluid obeying power law behavior is given by 

r|'=m(2n^y"''^^ (2.1) 

where m is the power law consistency index, n is the power law index 
and He is the second invariant of the rate of strain tensor, ey. 
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In order to complete the problem formulation, it is essential to have a 
mathematical model of cylinder-cylinder interactions. The spherical and 
cylindrical free cell models have been successful in predicting the macroscopic 
flow phenomena in particulate systems(Tripathi and Chhabra, 1 992). The two 
most extensively used models are the free surface model (Happel, 1958) and 
the Kuwabara model (1959). 

In this study, the Happel’s free surface cylindrical cell model is 
employed to study the flow behavior around a collection of cylinders. The free 
surface model consists of a cylinder of radius Tq surrounded by an imaginary 
cylindrical envelope of radius r* such that the porosity of each cell is equal to 
the porosity of the array of cylinders. Tlie cell boundary at r„ is frictionless, 
thereby implying the non-interacting nature of cells. In the Happel free surface 
model, a zero shear condition is imposed on the outer boundary. Depending 
on the distance r„, one can simulate different porosity conditions which 
defined by e = I - r„~^ . For the limiting case of a single cylinder in the fluid 
flow, porosity will become 1 . 


A cylinderical co-ordinate system is chosen for convenience. The angle 
0 is measured from the front stagnation point of the cylinder as shown in 
Figure 2.1(b), and To is the radius of inner cylinder and ro, that of the outer 
imaginary cylinder. The Fluid is approaching with a unifomi velocity of U and 
temperature . The components of the velocity in r- and 0- directions are 
and Vg respectively It is assumed that the flow is two dimensional and steady. 
The flow variables do not depend upon the coordinate z. 
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dimensional view 


2.2 Governing Equations for fluid flow 

With the above assumptions, the governing equations, namely the equations of 
continuity , motion and energy in cylindrical coordinates can be written as; 

The velocity components in the r, 6, z directions are 

Vr~Vr(r.e) ve'=ve(r,e) v,'=0 (2.2) 
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Similarly, the temperature can also be expressed as Tg(r',9). 

Under these conditions, the equation of continuity for an incompressible fluid 
is given as ; 


l acrv/) ^ 0 

r’ dr r' 39 


(2.3) 


and the components of the equations of motion are written as ; 
/•'-component : 


^ 1 rXrv/^) ^ 1 d{v\vg') 
dt' r' dr' r' 39 r' j 


dp' 

J? 


^ 1 d {r ^ 1 dr ^ £_oo_^ 
r' dr' r' 39 r' ) 


(2.4) 


0-component 


^'o , ^ ^(^0 ) , ^ ^(^'r ^ fl ) I '>^'r ^9 




dt' r' 39 


dr' 


r' 39 


+ 


1 dif'^^ro) I ^ ^’^'ee 


dr' r' 39 


(2.5) 


where x'y = 2 qVjj (i j = 0,r) 

where q' is the (non-Newtonian) viscosity 


( 2 . 6 ) 


q’ = nxny 

where O'e =8' Ve’^eo+s' V2(e' Ve’^9z+e'a^) 


In the above expressions 8ij is the rate of strain tensor which is related to 
the velocity components by 
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dr 


p. =1£^4.2 :c: 

r' dd r' 


(2.7) 

( 2 . 8 ) 


P* =: F* 

^ rQ ^ Or 


1 

r .^1 


1 dv/ 

— 

r' — 


+ ^ 

2 


U'j 

r' 36 _ 


(2.9) 


and for a power law fluid equation ( see 2. 1) is 

n' = m(2n’,;)^""*^^^ 


Now all the quantities in the above expressions are rewritten in their 
dimensionless form by using the following transformations. 


V,- 


U 


Vq 


u 



p = 


p~p' 

U2pU^ 



where T)^ = w 






( 2 . 10 ) 


Then, the continuity equation will remain unchanged and the equations from 
(2.4) and (2.5) remain unchanged in their form, and become 
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1 a(rv,) ^ I dvg ^ Q 
r dr r 39 


( 2 . 11 ) 


, 1 , 1 5(v,Va) Vfl^) 1 ap 

8 t r dr r 86 r ^ 2 dr 

^ 2” fl djrt^) ^ 1 dx,o Tqo^ 
Re(r dr r dO r ) 


(2.12) 


, 1 ^(H)') , 1 ^(^r^o) , v.Vq 
di r 86 r dr r 

^v{\ a(rx,) ^ iar„„^ 

Re(r’ dr r 89 j 

Substitution of the equations (2.6) through (2.9) into equations (2.12) 
and (2.13) yields 


1 1 dp 

iTm 


^ 1 g(rv/) ^ I 8{v,Vo) v/) I dp ^ l^pf d f I 5(rvj ) ^ 1 d\ 2 dvg '' 

^dt' r dr r 89 r J 2 dr Re (arfr dr J r^ 99^ r^ 89 ^ 

Re I " dr r 99, 


(2.14) 

dv„ 1 a(Vo^) 1 d(rv^Vg) \ \ dp ^ 2''p f d fl 5(rv0) ^ ^ 1 9\9 ^ 2 dv^ ^ 

^ dt ^ r dO ^ r dr ^ ^ 2 r d9 Rei^Sri^r dr j r^ d9^ r^ 99 ^ 

2‘"*'>r dp £<^,57}] 

Re [ ” ar r 99 , 

(2.15) 


where the Reynolds number, Re = 


pUM(2r^y 

m 


(2.16) 


Equation 2.12 and 2.13 can be solved numerically to obtain the velocity and 
pressure field. 


Furthermore, it is customary to introduce the corresponding vorticity function 



as; 


03 = 


dr 


r 86 r 


(2.17) 


the above equation become 


= + + — fa ^ + £^^1 

2 dr drf 2 J ** Re r 56 Re v ^ 5r r 86 j 


(2.18) 


15/? I5fv/+v''l l^rido) 2"*'( dri dr}'] 

2r dO r 56 2 J Re 5r Re v 5r r 50 j 


(2.19) 


These equations can be used to specify the boundary conditions for the 
pressure. 

2.3 Boundary Conditions 

In the Happel's free surface model, the fluid is assumed to be stationary 
and the cylinder as moving. The no-sIip boundary condition is employed on the 
cylinder surface(at r =1). For the cell surface(r = r„), the outer radial velocity 
and the shear stress is taken as zero. 

In Kuwabara's model, a zero vorticity is imposed on the outer surface 
instead of zero shear stress.Thus the boundary conditions are. 


atr-1 Vr=-cos0 V0=sin0 (2.20) 

at r=roo Vr=0 TrS^O in Happel Model, O)re-0 in Kuwabara Model. 

( 2 . 21 ) 

By a Gallilean transformation, we can convert this system to one where 
the cylinder is stationary and the fluid is moving, i.e., 

atr=/ v^=0 V0=O (2.22) 

at r=roo Vr=-cos0 t,. 0=O in Happel Model, <»r^0 in Kuwabara 
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Model. 


(2.23) 


These latter boundary conditions are used in the flow computations in this 
work. 

Fig 2.2 Cell model 


A. 


D 



2.3.1 Cell model and boundary conditions 


Boundary 

Boundary conditions 

AB 

vo=«f'=a|^=o 

ad 

BC 

11 

11 

<r> 

; > 

1 11 
! ^ 

CD 

dO 

AD 

Vr=-cosB, 'F=-rsin0, Xre =0 (Free 

surface cell model) and 

, (»^e=0(Kuwabara model) 


12 













Iln addition to the boundary conditions out lined in eq. (2.20) & (2.21), the 
condition of axisymmetry was used at 0 = 0 and B =n planes. 


2.4 Energy equation for Heat Transfer 

the equation of energy in terms of the transport properties, for non-Newtonian 
fluids of constant p,C^and k , and neglecting the viscous dissipation term, can 

be written as; 


pc. 


dP 


dirv; 1") 1 a( VoT') 

. ^ - 


dr' 


89 


K 


r' dr 


j df 
" dr' 


1 d^T' 

,'2 QQl 


(2.24) 


where p,Cp and k are the density, heat capacity and thermal conductivity of 
the fluid respectively which are assumed to be independent of temperature. 

Now all the quantities in the above expressions are rewritten in their 
dimensionless form by using the transformations as outlined in eq.(2.10) 

(not for temperature). 

In heat transfer applications, it is a common practice to employ two 
different type of thermal boundary conditions, namely, a constant temperature 
( j;) or a constant heat flux(9,, per unit length) at the surface of solid cylinder at 

r=l. 

Thus, depending upon the type of boundary condition, the non-dimensional 
temperature can be defined accordingly. Thus ,for the case of a constant 
temperature at r = 1, one can introduce the following dimensionless 
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temperature 7]: 



The non-dimensionalised energy equation does not depend upon the way the 
temperature is non-dimensionalised in these two cases. The only difference is 
that the boundary condition is changed at the cylinder surface i.e. at r = 1 . 

2.5 Boundary conditions for energy equation 
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As mentioned previously, two cases are considered here. For the case of 
uniform temperature at the solid cylinder : 


atr=l;T=l (2.32) 

atr=r„;T=0 (2.33) 

on the otherhand, when a constant heat flux is maintained at inner solid 
cylinder: 

at r= I; -^ = -1 (2.34) 

dr 

atr=r^;T=0 (2.35) 


Boundary 

Boundary condition 

AB 

^ = 0 
de 

BC 

dT 

— = -1 (for constant heat flux) 
dr 


T = -1 (for constant inner cylinder temperature) 

CD 

^ = 0 


dQ 

AD 

T = 0 


Table 2.2: Boundary conditions applied to figure 2.2 for energy equation 


2.6 Calculation of Nusselt number 

2.6.1 Constant surface temperature condion 

At the surface of the solid cylinder ( r = 1) equating the rate of heat transfer 


^suifact 


-K 


dP 

dr' 




Jr'*r0 


(2.36) 




Introducing dimensionless variables, 




dr 


HT,-To) 


(2.37) 


Nu^ 


h{2r,) _ , ^ 

^ L dr X., 


(2.38) 


2.6.2 Constant heat flux condition 

Again, one can equate the rate of heat transfer at the surface of the solid 
cylinder ( r = I ),as : 


ar' 

17 


-KT'lrur ~^o) (by defmation) 


Nov^ introducing dimentionless variables, 


(2.39) 


K Sr, ^,^1 


ro dr 


= hTl, 


K 


above eq.(2,39) gives 


dr 


r»i 


(2.40) 


(2.41) 


equating the first and third term of eq.(2.39) gives 


ro dr 


= hT 

r»l r»I 


(2.42) 


Putting the dimensionless temperature gradient from eq. (2.40), (2.41) gives, 


K 


■■hT 


r»l 


(2.43) 
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or 


(2.44) 


Nu„ 


H2r,) 

K 


rati 


The Nusselt number so evaluated relate to a fixed value of 9 ,and it is useful to 
Define an average Nusselt number Nu defined as follows: 


NUg 



(2.45) 


Thus, equations (2. 14), (2. 15) and (2.31) provide the theoretical 
framework for calculating the complete flow and temperature fields. These 
equations subjected to the boundary conditions have been solved numerically 
using the finite difference method.These results can be post processed to infer 
the values of the Nusselt number as a function of the Reynolds number and 
Prandtl number as well as porosity for the two types of thermal boundary 
conditions used here in. 
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Chapter 3 


NUMERICAL FORMULATION AND SOLUTION 

PROCEDURE 

The energy equation is solved using the finite difference technique.The 
equations governing fluid flow are solved using the finite difference method used 
in a previous work by Shibu (2000) and Shibu et al. (2001) and these solutions are 
used as input to the energy equation. The major difficulty encountered in the solving 
the Navier-Stokes equation, is due to the non-availability of any obvious equation 
which can be directly solved for the pressure. Therefore, the pressure field is found 
indirectly by finding a velocity field that satisfies continuity, along with the Navier- 
Stokes equations. 'Hiis involves iterations of the pressure and velocity fields. Such 
iterations have a very poor convergence on a normal grid. As a remedy, usually a so- 
called “staggered grid” is used for solving the equations. Therefore the energy equation 
is also solved using a staggered grid in this work. 

In the previous work, the SMAC-Implicit scheme was employed for solving the 
Navier-Stokes equations. This method is a pressure corrector method, which then uses 
a Poisson equation for the pressure corrections to obtain correction of the velocity 
fields, and the “pseudo-pressure” field which enforces the continuity equation. The 
pseudo-pressure, however, does not satisfy the correct pressure equation which has to 
be obtained by solving another Poisson equation with the correct boundary conditions. 
In this work, the implicit scheme with Gauss-Seidel iterations and under-relaxation is 
used for the solution of energy equation to obtain the temperature field. Although only 
the steady state solutions are required, we use transient equations to time-step from 
arbitrary initial conditions to the steady state. 
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3.1 Staggered grid 


In a staggered grid arrangement the computational domain is divided into 
number of cells as shown in Fig 3.1. The locations of scalars (pressure, viscosity and 
temperature) are at the centre of the cells, while those of vectors(velocities in r- and 9 - 
directions) are at the centre of cell faces to which they are normal. In such an 
arrangement, pressure difference between two adjacent cells directly affect the FD 
equations for the velocity components located at the interface of these cells. Thus the 
pressure and velocity corrections are directly related to each other, and the pressure- 
velocity iterations converge quickly. 

Cells are labeled with index (i,j) which denote the cell number as counted in 0 -and 
r-direclions, respectively; thus Ty is the temperature at the centre of cell (ij), while 
Vr(ij+l/2) is the radial velocity at the centre of the face between cell (ij) and cell 
(i,j+l), and so on. Similarly Vfl(i+1/2J) is the $ - component of velocity at the centre of 
face between cell (ij) and cell (i+l,j), and so on. Because of the staggered grid 
arrangements, the velocities are not defined at the cell centers, but whenever required, 
they are found by simple interpolations. Boundary conditions are imposed by setting 
appropriate velocities in the fictitious cells surrounding the physical domain (see 
Figure 3.2). In other words, the fictitious cell velocities are expressed in terms of 
boundary cell velocities for the implementation of boundary conditions. 

In Fig 3.2 B’C’F’G’ shows the actual physical domain. The domain ADEH is 
the computational domain. Fictitious cells are shown by dotted lines. Cells CDC’D’, 
EFE’F’, ABA’B’ and GHG’H’ are numbered as (1,1), (m,l), (l,n) and (m,n) 
respectively, where m and n are the number of cells in 9 -and r- directions respectively 
including the fictitious cells. 



I*ig 3. 1 'I'wo Dimensional Staggered Grid showing the locations of descretised 

variables. 



Fig 3.2 Grid Arrangement 




3.2 The Energy equation 


Consider the equation (2.3 1) governing the flow of heat for a power-law fluid 
Formulated in chapter 2, rewritten here 


dT ^ 1 din^J) ^ 1 d{v,T) 
dt r dr r dO 


RePr 


r dr 


f 


dr 


+ - 


1 

r^ d9^ 


(3.1) 


where Re, Pr are the Reynolds number and Prandtl number, respectively, and 
, Vg ,T stand for the non - dimensional velocities in r - and d - directions, and 
temperature, respectively. 


Eq. (3.1) can be slightly rearranged as : 


V dt r dr r dB j 


RePr 


vv 


aT 1 dT 

“ — i 

5/"* r dr 


\ 


+ 


1 d^r'' 
r^ d9\ 


(3.2) 


The above equation represents the equation of energy, in cylindrical co-ordinates. 
This equation consists of a transient term, convection terms and diffusion terms, whose 
discretisations are described in the ensuing sections. 


3,3 Discretisation of derivatives 


The center of discretisation for energy equation is (i, j). The discretisation of the 
derivatives of eq. (3.2) is carried out as follows, 

• Discretisation of the transient term 


1 

m#i+l *T*n 

l<v “ 

At 


• Discretisation of convection term containing 


1 d(rvj) 

_ 1 

r. 

j 



^ rpn-¥l , rpn-^l ^ 

-r jV J 

^2 ^2 

f' np-w+l , nr»rt+I 

> ■ 

[ 2 J 

1 2 J 

r dr 

Ar 
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Discretisation of convection term containing Vg 


r cO 


1 


V 

^i+\J 


^ ^TT/j+l , fpn+l 


^ I 

1 2 

[ z J 

Vi. 

1 2 J 



A0 


• Discretisation of second order derivative of temperature with respect to 
radius in diffusion term 


d'T 


jTfn + l “j“ 

df 

iJ 

j 

<3 

j 


* Discretisation of first order derivative of temperature with respect to 
radius in diffusion term 


I dT 

_ 1 

o-»n+J rpn-¥\ 

^ij+l ^ij-l 

r dr 

0 

2Ar 


• Discretisation of second order derivative of temperature with respect to 9 
in diffusion term 


1 d^T 

ft dQ^ 


ij 2 

0 


rpn^hl r\rpn-¥\ , n^AI + l 


3.3,1 Discretisation of Boundary conditions 

The fictitious cells shown in Figure (3.2) facilitate the discretisation of boundary 
conditions. 
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(a) For the constant cylinder surface temperature case 


Consider cells ... (i -1,2), (i,2), (i+1,2)... as cells on the inner boundary C’F’. So cells 
...(i-1,1 ), (i,l ), (i-‘ 1,1)... are fictitious cells along D’E’. So when we apply boundary 
condition for T at inner surface of the cylinder it is implemented as follows ; 

T'o = W^urfacc) - T,. (for i = 2 to m-1) 

Similarly, consider cells ... (i - 1, n - 1), (i, n-1), (i+1 ,n-l) ... as cells on the outer 
boundary' B'G’. So cells ... (i-l,n), (i, n), {i+l,n) ... are fictitious cells along A’H’. 

So when we apply the zero- derivative boundary condition for T at the outer 
boumiary, it is as follows : 

(for i=2to m-1 ) 

Along B'C' \ve have the discretisation of boundary condition as follows, 
similarly along E’G’ we have the discretisation of boundary condition as follows, 

(b) For the constant heat flux surface ( r = 1) 


Here only the boundary condition at r = 1 is different from the previous case. 

Consider cells ... (i-I, 2), (i, 2), (i-l, 2)... as cells on the inner boundary C’F’. So cells 
...(i-1, 1), (i, 1), (i+1, 1) ... are fictitious cells along D’E’. So when we apply boundary 

condition, — = -1 , for T at the surface of the solid cylinder as : 
dr 

^ , = (Ar) + (for i=2 to m-1) 

The boundary conditions elsewhere remain the same as in the case of the constant 
surface temperature. 
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3.4 Calculation of Nusselt number 

(a) For constant cylinder surface temperature 

The defining equation for Nusselt number for this case is given by Eq. (2.38) 
fonnulated in chapter 2. The local Nusselt number ( for a fixed value of 0 ) is calculated 
for each surface cell using the following grid point temperature, 

Ar 


(b) For constant heat flux at the surface (r = 1) 


The defining equation for Nusselt number for this case is given by Eq. (2.44) 
formulated in chapter 2. Here the Nusselt number is twice the inverse of the non- 
dimensionalised temperature at the surface. 


Nuii) 


I'he surface averaged values of the Nusselt number are obtained by simple 
Simpson’s one third rule. The values of Nusselt number are important output quantities 
which are plotted in the next chapter. 


3.5 Choice of Numerical Parameters 


The grid size used was typically 50x20 for systems with porosity ~0.4, 50x20 for 
porosity ~0.5, 60x40 for porosity ~0.6 and 60x100 for higher porosities (e >0.6). 

All results reported herein have been checked for mesh independence by using at 
least two different meshes. In broad terms, as the value of the power-law index 
decreased below the value of the unity, it became progressively more difficult to achieve 
the required level of convergence, especially for the values of n < 0.5. However in each 
case, the convergence criterion applied for the temperature field by Gauss-Siedel 
method was AT < 0.00001. And sometimes an under-relaxation factor of 0.7 was used 
in the Gauss-Siedel iterations to obtain the desired level of convergence. 

The CPU times were up to 100 hours on a pentinum 11 processor for the Navier- 
Stokes momentum equation solver, but just 1 hour for both constant temperature and 
constant heat flux cases for the energy equations solver.The time step used are in range 
of 0.01 - 0.0001 for the Prandtl number range of 1-500. 
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Chapter 4 

RESULTS AND DISCUSSION 


In this work, the numerical results of heat transfer for non-Newtonian fluids 
to tube bank systems have been obtained for the parameter ranges 0.40 < e < 0.6, 
1> n > 0.5, 1< Re < 500 and I < Pr < 500 for both boundary conditions, i.e., the 
constant wall temperature and constant heat flux at the surface of the solid 
cylinder. The . computer system used to perform these calculations are Linux 
systems of non-Newtonian fluid mechanics laboratory. Chemical Engineering 
Department, IIT Kanpur and servers in the Computer Centre working on Linux as 
their platform. This work mainly concentrated on high Reynolds number (Re > 
1.0) and at high Prandtl number (Pr^I.O), and the heat transfer rate in shear 
thinning fluids is seen to be always higher than that in Newtonian fluids. The 
results obtained in this work are new and as far as known to us have not been 
reported previously. 

4.1 Validation Of Newtonian Case Results 

Initially the solution procedures were validated by comparing the present 
results with the analytical results available in the literature for a Newtonian fluid. 
'I’hus the accuracy of the present results for Newtonian case can be established by 
comparing the present results with the published work of Mandhani et al. (2002). 
Mandhani ct al. (2002) have obtained the results of heat transfer, namely Nusselt 
number, for the flow of a Newtonian fluid for Pr=0.7 & Pr=7.7 over an array of 
cylinders. A comparison between these two sets of results is in Table 4.1. This 
table shows that the two values are close and the maximum divergence is less than 
3.5%. 
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l ahlc 4. 1 : C omparison of present work with the Mandhani et al. (2002) 


for constant surface temperature condition. 


€ 

\ 

Re 

Pr 

Nu(Present 

value) 

Nu ( Mandhani et al. (2002) 

0.99 

50.0 

0.7 

3.300502 

3.191410 

1 




7.7 

8.060382 

7.634112 

1 0.4 

500.0 

0.7 

22.468962 

22.601339 



7.7 

57.060432 

56.718174 

1 0.4 

1 

1 0.0 

77 

11.037856 

11.025297 


Nc.xt, tlie present solution procedure was validated by comparing the results 
for the case of single cylinder with the published work of Lange et al. (1998). 
Lange (1998) have obtained the results of heat transfer, namely Nusselt number, 
for the flow of air (Pr- 0.7148) over a single cylinder case. In the present case, the 
approximation to a single cylinder case was achieved by using a rather high value 
of porosity, namely 0.99, which makes the outer boundary ten radii away from the 
cylinder. The present and the published results are shown in the Table 4.2. The 
two values are seen to differ by about 10 % atmost. 

Table 4.2: Comparison of published work with present work for air 
( Pr - 0.7 1 48) for constant surface temperature condition 


Re 

Nu (Present value) 

Nu (Lange eta/. (1998) 

1 

0.913917 

0.816 

10 

1.665181 

1.8101 

20 

2.195308 

2.408678 

50 

3.296901 

3.642350 

100 

4.444632 

5.127775 

200 

1 6.416176 

7.420205 

500 

11.036654 

12.675766 
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Lange et fl/.(1998) correlated their numerical results by the following 
expression for the constant temperature boundary condition: 

Nu - OmiRe^-^ + 0.734 Re^ (4. 1 ) 

where x = 0.05 + 0.226 (4.2) 


These results apply for air in the range of 10''* < Re < 200 and entail an 
uncertainty of approximately 1.5 percent. The values presented by Lange et 
fl/.(1998) in their graphs and those calculated using Eq. (4.1) and (4.2) seem to be 
somewhat inconsistent. Limited results were also obtained for a single cylinder 
case by setting the porosity 6=0.999 which makes the outer boundary about 30 
cylinder radii large. For Re=100 and Pr=0.7, the present and literature values 
(Lange et a/. (1998)) are 4.42 and 4.5 respectively which differ by less than 2 % 
from each other. In this study, a much more fine grid was used than in the study of 
Lange et al. (1998). 

Furthermore, the present solution procedure was validated by comparing 
with an expression due to Leclair and Hamielee (1968). It was done by comparing 
the jii factor for porosity 0.99. The expression for factor is given in eq.(4.3) 
and eq.(4.4) respectively. The corresponding comparisons is shown in table 4.3. 




0.833 


1.05 Re 


- 0.582 



RePr*'^ 


(4.3) 

(4.4) 
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Tabic 4,3: Comparisons of j,i for porosity 0.99 for air Pr=0.7 


Rc 

j„ (Using eq.(4.3) 

(Using eq.(4.4) for 



present results) 

1 

l.(^88 

1.0206 

10 

0.2772 

0.1850 

.20 

o.issi 

0.1218 

j 50 

1 0.1086 

1 

0.0731 

1 100 ; 

! 0.0725 

1 

U— - - - - - 

0.0493 

^200 ' 

1 

1 0.0484 

0.0355 

' 500 i 

L -J 

0.0284 

0.0244 


The comparisons presented in Tables 4.2 and 4.3 illustrate the applicability 
of the present approach to an otherwise extremely difficult flow problem involving 
incompressible Newtonian fluids which represent a special case of power law fluid 
model (n=l). 

4.2 RESULTS FOR NON-NEWTONIAN LIQUIDS 

The new results for non-Newtonian fluids for porosities 0.4, 0.5 and 0.6 are 
given in 'fables 4.6 to 4.20. In general, it became increasingly difficult to obtain fully 
converged results at high values of voidage, lower values of n and / or higher values 
of the Reynolds number. 

4.2.1.(a) Results for Constant Surface Temperature 
The values of the mean Nusselt member for the constant temperature condition on 
the cylinder surface are presented in Tables 4.6 to 4.20 for a range of values of the 
Reynolds Number and voidage for 1 ^ Pr < 500 range and 0.5 < n ^ 1 range. For a 
constant value of power law index, as the value of porosity increases, the value of the 
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Niisselt number decreases for the same value of Reynolds number and Prandtl 
number. 

It is seen that the value of the Nusselt number increases as the power law 
index decreased, i.e., shear thinning lowers the resistance to heat flow and this 
facilitates heat transfer. This effect increases with the increasing values of the 
Reynolds number. Also, as expected the value of the Nusselt number is always larger 
for the case of constant heat flux condition than that for the constant temperature 
condition. The values of the local Nusselt number are also plotted as a function of 
Reynolds number and Prandtl number for different values of power law index in the 
fig. 4.1 to 4.52. Similarly the values of the local Nusselt number are plotted as a 
function of Reynolds number and Prandtl number with porosity as a parameter in the 
Fig. 4.103 to 4.136. 

An examination of these results suggest the following overall trends: 

(i) For a given porosity and values of n, the influence of Prandtl number 
increases with increasing values of the Reynolds number. This is 
presumably due to the growing thermal boundary layers at high 
Reynolds numbers. 

(ii) At low values of Reynolds number, the results are influenced much 
more by the porosity than that by the Prandtl number. However, as the 
values of the Reynolds number progressively rises, Prandtl number and 
shear thinning exert greater and greater influence in determining the 
value of Nusselt number. 

(iii) And also the results, at low values of Reynolds number, are not much 
influenced by power law index. However, as the value of the Reynolds 
number progressively rises, shearing action influences the heat transfer 
rate. 

The inspection of Fig. 4.1 to 4.52 and 4.103 to 4.136 suggests that at low Peclet 
number (= RePr), the Nusselt number almost has a constant value over the entire 
surface (Fig. 4.1). This is purely a conduction result. As the value of the Prandtl 
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number and/or Peclet number rises, the value of the Nusselt number is seen to be 
maximum at the front stagnation point, it decreases along the cylinder surface, 
attains a minimum near the rear stagnation point and finally shows an upturn in 
some cases (Fig. 4.2-4.50 and 4.103 to 4.136). As the value of the Reynolds number 
increases, the minima in the Nu-0 plots become more pronounced which is also 
seen to depend strongly on the value of the porosity. In dense systems (low values 
of porosity ), the Nusselt number shows larger variation over the surface. When the 
Reynolds number is greater than 100, the value of the Nusselt number seems to peak 
again after the minimum value. This is presumably due to the interference of the 
wakes in concentrated systems. Finally, the value of 9 at which the Nusselt number 
goes through a minimum value keeps moving forward which is in the line with the 
flow field. Once again we can observe that the Nusselt number increases with shear 
thinning action along the cylinder surface. The behavior of the Nusselt number 
variation along the cylinder surface for non-Newtonian fluids is qualitatively similar 
to that for a Newtonian fluid over an array of circular cylinders. 


4.2.L(b) Temperature Field 

Some further insights into the nature of heat transfer to/ffom non- 
Newtonian fluids can be gained by examining the isotherm plots. These are 
presented in figures 4.158-4.397, for constant temperature condition and the flow 
is from left to right in these figures. The inner circle shows the surface of the solid 
cylinder (r =1). These curves are plotted for temperature values of 0.01 to 1.0 in 
step of 0.01, the temperature being unity at r=l and zero at r= r„ . 
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It can be seen from fig 4.158-4.240., that for low values of Reynolds 
number and Prandtl number isotherms show front and rear symmetry for a given 
value of porosity and power-law index. This is due to the fact that heat transfer 
occurs mainly by conduction with negligible convection. As the value of Reynolds 
number and/or Prandtl number (i.e., Peclet number) keeps on increasing, this 
symmetry gradually disappears indicating the increasing contribution of 
convection as shown in figures 4.158-4.240 Same phenomenon is repeated with 
different values of porosity as shown in figures 4.241-4.397. In general, higher is 
the porosity, lower is the value of the Reynolds number up to which the 
temperature field displays front and rear symmetry. 

4,2.2.(a) Results for constant heat flux condition 


'I’he values of the mean Nusselt number for constant flux boundary 
condition are presented in Tables 4.6 to 4.20 for a range of values of the system 
parameters. As expected, the values of the mean Nusselt number are always higher 
under this condition than those for constant surface temperature condition, though 
the difference between the two values increases with decreasing porosity and/or 
increasing Reynolds number. It is seen that the Nusselt number increase as the 
power law index decreased, i.e., shear thinning lowers the resistance to heat flow. 
As the value of porosity becomes higher the value of the Nusselt number 
decreases for the same value of Reynolds number and Prandtl number, while for a 
given porosity the value of the Nusselt number increased both with Reynolds 
number and Prandtl number. In this case the local as well as mean Nusselt number 
are consistently higher than for constant surface temperature condition. The values 
of the local Nusselt number are also plotted as a function of Reynolds number and 
Prandtl number with porosity as a parameter in the figures 4.52-4.102 and 4.132- 
4.155. 
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1 he trend in the value of Nusselt number is same as that discussed for 
constant surface temperature condition. Here the value of the Nusselt number is 
higher at the front and rear stagnation point as are the mean values. Once again we 
can observe that the Nusselt number increases with shear thinning action along the 
cylinder surface. The behavior of the Nusselt number variation along the cylinder 
surface for non-Newtonian fluids is qualitatively similar to that for a Newtonian 
fluid over an array of circular cylinders. 

4.2.2.(b) TEMPERATURE FIELD 

Isotherm plots 

Some further insights into the nature of heat transfer of non- 
Newtonian fluids can be gained by examining the isotherm plots. These are 
presented in figures 4.398-4.473 for constant heat flux condition and the flow is 
from left to right in these figures. The inner circle shows the surface of the solid 
cylinder. These curves are plotted for temperature values of 0.01 to 1.0 in step of 
0.05, the temperature being unity at r=l and zero at r= r„ . 

It can be seen from fig 4.398-4.473 that for low values of Reynolds number 
and Prandtl number isotherms show front and rear symmetry for a given value of 
porosity. This is due to the fact that heat transfer occurs mainly by conduction 
with negligible convection. As the value of the Reynolds number and/or Prandtl 
number (i.e., Peclet number) keeps on increasing, this symmetry gradually 
disappears indicating the increasing contribution of convection as shown in figures 
4.398-4.631 As Peclet number (RePr) keeps on increasing it results in less number 
of isotherms for a definite value of porosity. The outer most isotherm is having a 
value of 0.05. Inspection of Figures 4.473-4.631 reveals that as porosity increases, 
Peclet number being the same, more and more isotherms are present as the system 
goes towards the sparse systems. Same phenomena can be seen for a different 
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value of porosity as shown in figures 4.473-4.631. In general, higher is the 
porosity, lower is the value of the Reynolds up to which the temperature field 
displays front and rear symmetry. 

4.3 Validation of non-Newtonian fluids results 

Initially, the present solution procedure was validated by comparing j fj 
factor for the special case of porosities of 0.75 (n=0.748) and 0.714 (n=0.83) with 
the experimental results of Adams and Bell (1968). The values presented by 
Adams and Bell ( 1 968) in graphs and those calculated using the equation: 


,/„ =0.599 Re ■‘’■“V 0.003 ( 4 . 5 ) 

for e =0.714, n=0.83 and Pr=740.56 (0.5%CMC solution). From the 
geometric details of tube bundles used by Adams and Bell (1968), their results 
have been recalculated in the form required here (as shown in Appendix A). A 
comparison of these results with the present numerical results is shown in the 
Table 4.4. 


Table 4.4. Comparison of the J // factors 


Re 

Nu 


J H numtricQl 

J H exp tl 

Jm 

numerical 

Jh 

-'^exptl 

5.0 

24.887276 

16.724 

0.122426 

0.091807 

1.3335 

7.5 

1 29.345327 

25.086 

1 0.096240 

0.070124 

1.37243 

10.0 

33.284916 

33.448 

0.081868 

0.057933 

1.42567 

20.0 

52.959656 

66.896 

0.0651305 

0.036598 

1.779619 
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Furthermore, u=0.75, n=0.748 and Pr=92L36 (1%CMC solution). From the 
geometric details of tube bundles used by Adams and Bell (1968), their results 
have been recalculated in the form required here (as shown in Appendix B). A 
comparison of these results with the present numerical results is as follows; 


Table 4.5. Comparison of the j factors 



Nu 




^ numerical 

Re 


J H numerical 

J H exp il 

Jh 

1.0 

13.823075 

3.7860 

0.4518098 

0.306785 

1.47272 

2.5 

19.361176 

9.4654 

0.2531237 

0.162835 

1.55448 


'“15442535"" 

18.9308 

0.1663151 

0.099594 

1.66793 


30.383961 

! 28.3962 

0.1324111 

0.073498 

1.80156 

'ibTo 

44.131569 

37.8616 

0.104240 

0.059164 

1.76186 


An examination of this Tables 4.4 & 4.5 shows that the two results differ by 
a maximum of 78%. While this may seem like a rather large discrepancy, but it 
must be borne in mind that finite wall effects and maldistribution of flow are 
always present in experimental study which are altogether neglected in the 
numerical simulations. Also, while the numerical simulations are based on 
constant thcnno-physical properties, significant temperature intervals are present 
in the study of Adams and Bell (1968). Bearing in mind these factors, the 
correspondence seen in Tables 4.4 & 4.5 are acceptable and satisfactory. 
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Table No.4.6; Reynolds number =1.0, porosity = 0.4 


11 

Pr 

Nu(t) 

Nu(f) 

I.O 

1 

1 

7.934439 

7.936596 

1 fo 

8.036373 

8.171126 

1 50 

9.545625 

10.459119 

100 

12.021615 

13.234608 

500 

21.006881 

21.82945 




0.8 

0.6 

1 

7.934540 

7.935972 

10 

8.036493 

8.166099 

50 

9.633212 

10.522439 

100 

12.093627 

13.262397 

500 

21,021563 

22.465096 




1 

7.960410 

7.961659 

10 

8.055016 

8.170110 

50 

9.780549 

10.614179 

100 

12.106388 

13.308061 

500 

21.147916 

23.050785 




0.5 

1 

7.960517 

7.961791 

10 

8.055350 

8.174158 

50 

9.784709 

10.620215 

100 

12.166709 

13.461456 

o 

o 

21.342908 

23.923456 
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I able No.4.7; Reynolds number =10.0, porosity = 0.4 


!1 

Pr 

Nu(t) 

Nu(f) 

I.O 

1 

8.0163260 

8.070533 

10 

12.143661 

13.316152 

50 

21.131322 

22.901655 

100 

27.0191656 

29.1559647 

500 

48.896416 

53.531244 




0.8 

1 

8.0307695 

8.155944 

10 

12.161279 

13.423942 

50 

21.244217 

23.234256 

100 

27.152752 

29.737281 

500 

49.281631 

53.725784 




0.6 

1 

8.035966 

8.162755 

10 

12.186206 

13.444005 

50 

21.340579 

23.3405218 

100 

27.2421238 

29.821318 

500 

49.760937 

53.827808 




0.5 

1 

8.059018 

8.186030 

10 

12.189671 

13.447103 

50 

21.396200 

23.396727 

100 

27.405266 

29.825638 

500 

49.84706 

53.959251 
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Table Ni).4.8; Reynolds number = 100.0, porosity = 0.4 


II 

Fr 

Nu(t) 

Nu(f) 


I 

12.900322 

14.420811 


10 

31.127928 

31938839 

1.0 

50 

59.168449 

63.912827 


100 

122.544151 

128.039513 


500 

213.221893 

218.962402 






1 

11324706 

14.942348 


10 

33.801826 

37.022709 

0.8 

50 

68.053444 

67.731087 


100 

124.526047 

132.804657 


500 

218.111557 

227.731079 






f 

14415411 

15.958244 ‘ 


^ 10 

38.623264 

42.873165 

0.6 


66T444885“ 

78.850395 


100 

135.347931 

157.943436 


500 

246.147354 

249.481934 






1 

14.754089 

16.772413 


10 

42.743153 

47.799911 

0.5 

50 

89.462036 

101.665321 


100 

i 45.083099 

168.832169 


500 

257.806702 

268.358063 
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'Fable No.4.9; Reynolds number =200.0, porosity = 0.4 


11 

Pr 

Nu(t) 

Nu(f) 

1.0 

n 

17.715052 

19.712749 

iO 

45.128601 

49.185493 

50 

113.196922 

122.697411 

100 

155.446106 

173.965134 

O 

O 

253.412811 

259.700409 




0.8 

1 

19.072891 

21.266714 

10 

50.811386 

55.656525 

50 

12~6. 156776 

147.413635 

100 

175.413345 

192.377274 

500 

267.338562 

276.814148 




0.6 

1 

20.898769 

23.610092 

10 

62.935993 

75.452873 

50 

153.894180 

171.878906 

100 

207.5146 

221.171646 

500 

278.584595 

304.141876 




0.5 

i 

22.320176 

25.352266 

10 

68.229858 

77.775620 

50 

169.132690 

187.320923 

100 

223.761124 

238.929703 


- 

- 

HHIH 
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I able No. 4. 10; Reynolds number =500.0, porosity = 0.4 


n ■ 

Pr 

Nu(t) 

Nu(f) 

1 

I 

29.364960 

32.431778 

1.0 

10 

71.549950 

80.174782 

50 

199.541245 

211.599335 

100 

242.640045 

254.563065 

500 






0.8 

1 

32.872173 

37.430359 

10 

95.310158 

114.532 

50 

227.004105 

239.641418 

100 

238.460571 

278.567413 

500 






0.6 

1 

37.485847 

39.1304206 

10 

432.6194154 

162.072034 

50 

235.441217 

247.321041 

100 

266.582794 

279.334762 

500 






0.5 

1 

38.006926 

41.832642 

10 

176.515778 

187.385422 

50 

256.983276 

268.021271 


272.207733 

279.483551 
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Table No.4. 1 1 ; Reynolds number =1.0, porosity = 0.5 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

1 

5.847853 

5.8487290 

10 

6.001798 

6.013905 

50 

7.985848 

8.6935014 

100 

10.201760 

10.991025 

500 

17.851946 

18.912055 




0.8 

1 

5.847913 

5.849944 

10 

5.993634 

6.054933 

50 

7.9951415 

8.719989 

100 

10.018479 

10.981187 

500 

17.417665 

18.892799 




0.6 

1 

5.848024 

5.850346 

10 

5.994559 

6.150870 

50 

7.96923488 

8.780470 

100 

10.1963233 

10.998605 

500 

17.457107 

18.693872 




0.5 

1 

5.8488031 

5.8610275 

10 

6.010251 

6.1725814 

EHHIli 

7.9964302 

8.783092 

DSEHIH 

10.2039752 

10.99879102 
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Table No.4.12; Reynolds number =100.0, porosity = 0.5 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

1 

j 

! 

! 

1 

10.877949 

12.233936 

10 

26.076292 

28.817455 

50 

48.445957 

53.074196 

Too 

64.332642 

69.672241 

500 

171.454330 

196.622559 




0.8 

1 

1 

f 1 

1 

(16 

1 

11.641509 

13.158134 

10 

29.342411 

32.247120 

50 

57.242017 

62.4138501 

Too 

71 . "549431 

78.942345 

500 

190.817993 

208.283539 




1 

11.90^25 

13.428152 

' 'iff 

31.567934 

34.684757 

"so 

60.582745 

66.015411 

'Too 

82.640030 

90.591469 

'506 

207.671860 

221.878677 




0.5 

1 

12.507766 

14.140772 

10 

34.966072 

38.579731 

50 

67.909447 

74.646080 

100 

109.641670 

133.203979 

500 

229.028152 

242.216324 
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I'ahle No.4. 1 3; Reynolds number =10.0, porosity = 0.5 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

i 

5.9967028 

6.16780162 

10 

10.00723109 

11.00173266 

50 

17.01846023 

18.97831062 

100 

21.89105602 

24.0965238 

500 

39.3967098 

42.4890116 




0.8 

1 

5.99^009 

6.168808 

10 

l”0.013325 

11.089046 

150 

17.353355 

19.082125 

100 

22.091156 

24.174515 

500 

'39.i24084 

42.592155 




0.6 

1 

1999010 

6.168867 

1 io 

10.023359 

11.099346 

50 

17.391918 

19.117094 

100 

22.143450 

24.218283 

!500 

39.481960 

42.625042 

1 



0.5 i 

1 

5.999635 

6.169416 

10 

10.0694650 

11.154601 

50 

17.545757 

19.293976 

100 

22.366032 

24.473021 

500 

39.895378 

42.993656 
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Table No.4. 14; Reynolds number =200.0, porosity = 0.5 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

1 

14.86104452 

16.98707740 

[o 

40.09704638 

43.44088601 

^0 

7^8826041 

78.9202765 

100 

122.72681053 

137.8510472 

500 

231.9778825 

242.9820344 




0.8 

I 

16.154152 

18.023235 

To 

42.264950 

46.192909 

50 

77.139107 

85.602890 

Too 

13^725281 

171.480865 

500 

256.684753 

265.087982 




0.6 

^ 1 

17.766781 

19.927652 

To 

50.340824 

55.811890 

i50 

91.730209 

105.045273 

100 

175.410706 

191.164352 

500 

291.873291 

293.46003 




0.5 

1 

18.866871 

21.247288 

10 

56.668926 

63.287262 

50 

138.407745 

159.342453 

100 

192.104858 

208.570724 

500 

293.558350 

304.239502 
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I able No. 4. 1 5; Reynolds number =500.0, porosity = 0.5 


!1 

Pr 

I 

Nu(t) 

Nu(f) 

! 

1 

1.0 

1 . 

24.671692 

27.197191 

10 

61.873554 

67.928024 

50 

164.440140 

187.232742 

100 

234.197708 

248.330505 

500 






! 

0.8 

1 

27.'427446 

30.486328 

10 

78.516046 

132.779495 

50 

20^495255 

213.182373 

100 

248.166275 

259.812378 

500 






0.6 

1 

30.782507 

34.594124 

10 

99.882957 

117.319771 

50 

227.452621 

242.090790 

100 

283.049622 

293.547577 

500 






0.5 

1 

32.213848 

36.280319 

lEHi 


145.886871 

50 

237.430145 

250.633865 

DSHBII 

291.745544 

298.755615 

E^HHI 
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Tabic No.4.16; Reynolds number =1.0, porosity = 0.6 


n 

Pr 

Nu(t) 

Nu(f) 

i.O 

1 

4.414944 

4.417809 

10 

4.623365 

4.809139 

50 

6.647354 

7.25817 

100 

8.348444 

9.108577 

'500 

14.361089 

15.454818 




0.8 

1 

4.414996 

4.417840 

10 

4.626966 

4.819084 

50 

6.651194 

7.295737 

100 

8.342834 

9.123822 

^00 

14.418490 

15.669935 




0.6 

0.5 

1 

4.424533 

4.427312 

10 

4.633632 

4.823524 

50 

6.639102 

7.275495 

100 

8.376125 

9.1339940 


14.431066 

15.746615 




1 

4.424587 

4.427317 

IQflUII 

4.634986 

4.829026 

50 

6.641713 

7.277069 

100 

8.388171 

9.208723 

500 

14.536267 

15.941373 
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I'ahlc No.4, 1 7; Reynolds number =10.0, porosity = 0.6 


11 

Pr 

Nu(t) 

Nu(f) 

I.O 

0.8 

1 

4.633389 

4.839238 

10 

8.448704 

9-380105 

50 

14.611475 

16.155329 

O 

o 

18.572514 

20.453854 

500 

31872940 

35.723591 




I 

4.637657 

4.846606 

10 

8.451797 

9.388928 

50 

14.593246 

16.134794 

100 

l[8.y40773 

20.415014 

500 

32.828190 

35.750427 




0.6 

0.5 

1 

4.638368 

4.847233 

10 

8.484939 

9.399450 

50 

14.490101 

15.896059 

100 

18.585667 

20.478066 


32.838493 

35.816815 




1 

4.638368 

4.8510602 

10 

8.5010790 

9.41087662 


14.6125062 

15.9055468 

lESIHHi 

18.610002317 

20.50011567 

BSBHi 

32.8409201 

35.8904272 

HHHI! 
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Table No.4.18; Reynolds number =100.0, porosity = 0.6 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

1 

9.383591 

10.586508 

10 

22.463772 

24.735071 

50 

41.510197 

45.189659 

100 

54.223763 

58.949402 

500 

138.227753 

158.977188 




0.8 

1 

9.714847 

10.736677 

10 

24.108482 

26.395256 

50 

45.231358 

48.956703 

100 

59.765381 

64.620941 

500 

159.776733 

169.460297 




0.6 

1 

9.919545 

11.222487 

10 

26.787050 

29.408550 

50 

51.976742 

56.860264 

100 

71.358162 

79.854462 

500 

174.877899 

189.156952 




0.5 

1 

10.099916 

11.416263 

10 

26.788191 

29.482941 

50 

52.206837 

57.288868 

100 

72.434319 

82.889786 

500 

176.594315 

192.441620 
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Table No.4.19; Reynolds number =200.0, porosity = 0.6 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

1 

13.108223 

14.641165 

10 

32.736610 

35.696430 

50 

61.060459 

66.398125 

100 

101.986938 

126.515434 

500 

208.733551 

223.521561 




0.8 

1 

14.027789 

15.648171 

10 

36.270531 

39.614288 


70.541786 

77.570686 

100 

129.566742 

137.270782 


236.029968 

249.409546 




0.6 

1 

15.921167 

17.148172 


42.826439 

45.091190 


85.816048 

91.924561 

100 

119.369812 

134.720230 


253.720093 

265.709259 




0.5 

1 

16.435616 

18.387285 


47.680515 

52.78582 


87.990746 

94.542900 


162.498703 

179.779251 

500 

278.983856 

286.505707 
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Table No.4.20; Reynolds number =500.0, porosity = 0.6 


n 

Pr 

Nu(t) 

Nu(f) 

1.0 

1 

22.080944 

23.800337 


56.175209 

60.045734 

BHHIi 

137.498016 

172.751877 

100 

188.836365 

219.978104 


278.377197 

287.058685 




o 

bo 

1 

24.010752 

26.607241 

10 

66.757576 

75.825409 

50 

166.523087 

182.637512 

100 

221.351349 

233.809479 

500 

278.377319 

312.638275 




0.6 

1 

26.816025 

29.964300 

10 

80.452332 

96.514900 

50 

199.987335 

211.520767 

100 

257.210480 

260.914978 

500 

278.377441 

334.940674 




0.5 

1 

27.913700 

31.291368 

10 

81.723160 

91.012787 

50 

208.022507 

220.703491 

100 

268.975189 

271.482666 
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Angle in degrees 

Figure 4.1: Variation of Nusselt number with angle for Re=1.0,e=0.4 and n=1.0 
for constant surface temperature condition 



Angie in degrees 


Figure 4.2: Variation of Nusselt number with angle for Re=10.0,e=0.4 and n=1.0 
for constant surface temperature condition 
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Figure 4.3: Variation of Nusselt number with angle for Re=:l()0.0,e=().4 and n=1.0 
for constant surface temperature condition 
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Figure 4.4: Variation of Nusselt number with angle for Re=200.0,e=0.4 and n=1.0 
for constant surface temperature condition 
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Figure 4.5: Variation of Nusselt number with angle for Re=500.0,e=0.4 and n=1.0 
for constant surface temperature condition 
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Figure 4.6: Variation of Nusselt number with angle for Re=1.0,e=0.4 and n=0.8 
for constant surface temperature condition 
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Figure 4.7: Variation of Nusselt number with angle for Re=10.0,e=0.4 and n=0.8 
for constant surface temperature condition 
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Figure 4.8: Variation of Nusselt number with angle for Re=100.0,e=0.4 and n=0.8 
for constant surface temperature condition 
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Figure 4.9: Variation of Nusselt number with angle for Re=200.0,e=0.4 and n=0.8 
for constant surface temperature condition 
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Figure 4.10: Variation of Nusselt number with angle for Re=500.0,e=0.4 and n=0.8 

for constant surface temperature condition 
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Figure 4.11: Variation of Nusselt number with angle for Re=l(). 05 e =().4 and n=0.6 
for constant surface temperature condition 
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Figure 4.12: Variation of Nusselt number with angle for Re=100.0,e=0.4 and n=0.6 
for constant surface temperature condition 
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Figure 4,13: Variation of Niisselt number with angle for Re=200.0,e=0.4 and n=0.6 

for constant surface temperature condition 
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Figure 4.14: Variation of Nusselt number with angle for Re=500.0,e=0.4 and n=0.6 

for constant surface temperature condition 





Re=1, e==0.4,n=0.5(const wall temp) 
free surface model 


Pr=500.0 


Pr=100.0 


Pr=50.0 


Pr=10.0 


Angle in degrees 


Figure 4.15: Variation of Nusselt number with angle for Re=1.0,e=0.4 and n=0.5 
for constant surface temperature condition 
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Figure 4.16: Variation of Nusselt number with angle for Re=10.0,e=0.4 and n=0.5 
for constant surface temperature condition 
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Figure 4.17: Variation of Nvisselt number with angle for Re=100.0,e=0.4 and n=0.5 
for constant surface temperature condition 
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Figure 4.18: Variation of Nusselt number with angle for Re=200.0,e=0.4 and n=0.5 
for constant surface temperature condition 
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Figure 4.19: Variation of Nusselt number with angle for Re=L0,e=0.4 and n=(l5 
for constant surface temperature condition 
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Figure 4.20: Variation of Nusselt number with angle for Re— l,0,e=0,5 and n=0.8 
for constant surface temperature condition 
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Figure 4.23: Variation of Nusselt number with angle for Re=200.0,e=0.5 and n=0.8 
for constant surface temperature condition 
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Figure 4.24: Variation of Nusselt number with angle for Re=500.0,e=0.5 and n=0.8 
for constant surface temperature condition 
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Figure 4.27: Variation of Nusselt number with angle for Re=100.0,e=0.5 and n=0.6 
for constant surface temperature condition 
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Figure 4.28: Variation of Nusselt number with angle for Re=200.0,e=0.5 and n=0.6 
for constant surface temperature condition 
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Figure 4.29: Variation of Nusselt number with angle for Re=500.0,e=0.5 and n=0.6 

for constant surface temperature condition 
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Figure 4.36: Variation of Nusselt number with angle for Re=10.0,e=0.6 and n=1.0 
for constant surface temperature condition 


67 








Pr=500.0 


Re=100, 0=0.6, n=1-0(const waff temp) 
free surface model 


Pr= 100.0 


Pr=50.0 


Pr=10.0 


Pr=1 .0 


Angle In degrees 

Figure 4.37: Variation of Nusselt number with angle for Re=100.0,e=0.6 and n=1.0 

for constant surface temperature condition 
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Figure 4.38: Variation of Nusselt number with angle for Re=200.0,e=0.6 and n=1.0 

for constant surface temperature condition 
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Figure 4.40: Variation of Nusselt number with angle for Re=10.0,e=0.6 and n=0.8 
for constant surface temperature condition 
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Figure 4.41: Variation of Nusselt number with angle for Re=200.0,e=0.6 and n=0.8 

for constant surface temperature condition 
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Figure 4.42: Variation of Nusselt number with angle for Re=500.0,e=0.6 and n=0.8 

for constant surface temperature condition 
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Figure 4.43: Variation of Nusselt number with angle for Re=1.0,e=0.6 and n=0.6 
for constant surface temperature condition 



Figure 4.44: Variation of Nusselt number with angle for Re=10.0,e=0.6 and n=0.6 
for constant surface temperature condition 
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Figure 4.45: Variation of Nusselt number with angle for Re=100.0,e=0.6 and n=0.6 

for constant surface temperature condition 
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Figure 4.46: Variation of Nusselt number with angle for Re=200.0,e=0.6 and n=0.6 

for constant surface temperature condition 
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Figure 4.48; Variation of Nusselt number with angle for Re=1.0,e=0.6 and n=0.5 
for constant surface temperature condition 
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Figure 4.51; Variation of Nusselt number with angle for Re=200.0,e=0.6 and n=0.5 
for constant surface temperature condition 
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Figure 4.52: Variation of Nusselt number with angle for Re== 500 . 05 e=== 0.6 and n=0.5 
for constant surface temperature condition 
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Figure 4.55: Variation of Nusselt number with angle for Re= 100.0,6=0.4 and n=1.0 

for constant heat flux condition 
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Figure 4.56: Variation of Nusselt number with angle for Re=200.0,e=0.4 and n=1.0 

for constant heat flux condition 
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Figure 4.58: Variation of Nusselt number with angle for Re=1.0,e=0.4 and n=0.8 

for constant heat flux condition 
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Figure 4.63: Variation of Nusselt number with angle for Re=10.0,e=0.4 and n=0.6 
for constant heat flux condition 
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Figure 4.64; Variation of Nusselt number with angle for Re= 100.0, e=0.4 and n=0.6 

for constant heat flux condition 
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Figure 4.65: Variation of Nusselt number with angle for Re=200.0,e=0.4 and n=0.6 

for constant heat flux condition 


Figure 4.66: Variation of Nusselt number with angle for Re=500.0,e=0.4 and n=0.6 

for constant heat flux condition 
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Figure 4.67: Variation of Nusselt number with angle for R€=1.0,e=0.4 and n=0.5 ^ure 4.69: Variation of Nusselt number with angle for Re=100.0,e=0.4 and n=0.5 for constant heat flux concli 

for constant heat flux condition ! 
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re 4.70: Variation of Nusselt number with angle for Re=200.0,e=0.4 and n=0.5 for constant heat flux condi 

r 



Figure 4.68: Variation of Nusselt number with angle for Re=10.0,e=0.4 and n=0.5 

for constant heat flux condition 


84 




Pr=50.0 


Pr=10,0 


Re=500', e=:0'.4,n=0.5(const heat flux) 
free surface model 


Pr=1.0 


100 120 140 160 

Angle In degrees 


I 


Figure 4.71: Variation of Nusselt number with angle for Re=1.0,e=0 4 and n=0 5 

for constant heat flux condition 
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Figure 4.72: Variation of Nusselt number with angle for Re=1.0,e=0.5 and n=0.8 

for constant heat flux condition 
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Figure 4.73: Variation of Nusselt number with angle for Re=10.0,e=0.5 and n=0.8 

for constant heat flux condition 
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Figure 4.74: Variation of Nusselt number with angle for Re=100.0,e=0.5 and n=0.8 

for constant heat flux condition 
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Figure 4.81: Variation of Nusselt number with angle for Re=500.0,e=0.5 and n=0.6 

for constant heat flux condition 
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Figure 4.82: Variation of Nusselt number with angle for Re=1.0,e=0.5 and n=0.5 

for constant heat flux condition 
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Figure 4.85: Variation of Nusselt number with angle for Re=200.0,e=0.5 and n=0.5 

for constant heat flux condition 
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Figure 4.86: Variation of Nusselt number with angle for Re=500.0,e=0.5 and n=0.5 

for constant heat flux condition 
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Figure 4.87: Variation of Nusselt number with angle for Re=1.0,e=0.6 and n=1.0 

for constant heat flux condition 



Angle In degrees 


Figure 4.88. Variation of Nusselt number with angle for Re=10.0,e=0.6 and n=1.0 

for constant heat flux condition 
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Figure 4.91: Variation of Nusselt number with angle for Re=:500.0,e=0*6 and 0=1.0 

for constant heat flux condition 
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Figure 4.92: Variation of Nusselt number with angle for Re=10.0,e=0.6 and n=0.8 

for constant heat flux condition 
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Figure 4.95; Variation of Nusselt number with angle for Re= 1.0, 6=0.6 and n5=0.6 

for constant heat flux condition 



Figure 4.96; Variation of Nusselt number with angle for Ile=10.0,e=0.6 and n=0.6 

for constant heat flux condition 
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Figure 4.99: Variation of Nusselt number with angle for Re=500.0,e=0.6 and n=0.6 

for constant heat flux condition 
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Figure 4.100: Variation of Nusselt number with angle for Re=l,0,e=0.6 and n=0.5 

for constant heat flux condition 
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Figure 4.101: Variation of Nusselt number with angle for Re=10.0,e=0.5 and n=0.5 

for constant heat flux condition 
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Figure 4.102: Variation of Nusselt number with angle for Re=100.0,e=0.6 and n=0.5 

for constant heat flux condition 
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Figure 4.103: Variation of Nusselt number with angle for Re=200. 

for constant heat flux condition 
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Figure 4.104: Variation of Nusselt number with angle for Re=500.0,e=0.6 and n=0.5 

for constant heat flux condition 
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Figure 4.103:Variation of Nusselt number with angle for Re=l,n=0.8 and Pr=l 

for constant wall teperature 
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Figure 4. 104: Variation of Nusselt number with angle for Re=l,n=0.8 and Pr=10 

for constant wall teperature 
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Figure 4.106: Variation of Nusselt number with angle for Re=l,ns=0.8 and Pr=500 

for constant wall teperature 
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Figure 4.107:Variation of Nusselt number with angle for Re=10,n=0.8 and Pr=l 

for constant wall teperature 



Figure 4.108:Variation of Nusselt number with angle for Re=100,n=0.8 and Pr=l 

for constant wall teperature 
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Figure 4.109:Variation of Nusselt number with angle for Re=100,ir— 0.8 and Pr=10 

for constant wall teperature 
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Figure 4. 110: Variation of Nusselt number with angle for Re=100,n==0.8 and Pr=100 

for constant wall teperature 
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Figure 4.113;Variation of Nusselt number with anghf for Re=500,n=0.8 and Pr=50 

for constant wall teperaturc 
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Figure 4.114;Variation of Nusselt number with r nglc for Re-10, n=0.6 and Pr=l 

for constant wall tep 3 Jature 
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Figure 4.117:Variation of Nusselt number with angle for He“i00,n=0.6 and Pr=l 

for constant wall teperatiirt' 
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Figure 4.118:Variation of Nusselt number with angle for Re=100,n=0.6 and Pr=10 

for constant wall teperature 
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Figure 4.124:Variation of Nusselt number with angle for Re=l,n=0.5 and Pr=10 

for constant wall teperature 
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Figure 4. 126: Variation of Nusselt number with angle for Re=l,n=0.5 and Pr-100 

for constant wall teperature 
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Figure 4.127:Variation of Nusselt number with angle for Re=l,n=0.5 and Pr=500 

for constant wall teperature 



Figure 4.128:Variation of Nusselt number with angle for Re=100,n=0.5 and Pr=l 

for constant wall teperature 
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Figure 4.129:Variation of Nusselt number with angle for R(!=100,n=0.5 and Pr=10 

for constant wall teperature 
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Figure 4. 130: Variation of Nusselt number with angle for Re=100,n=0.5 and Pr=100 

for constant wall teperature 
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Figure 4. 131: Variation of Nusselt number with angle for Re=100,n=0.5 and Pr=500 

for constant wall teperature 
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Figure 4.132:Variation of Nusselt number with angle for Re=500,n=0.5 and Pr=l 

for constant wall teperature 
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Figure 4. 135: Variation of Nusselt number with angle for Re=l,n=0.8 and Pr=100 

for constant heat flux 



Figure 4. 136: Variation of Nusselt number with angle for Re=l,n=0.8 and Pr=500 

for constant heat flux 
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Figure 4.137:Variation of Nusselt number with angle for Re=10,n=0.8 and Pr=l 

for constant heat flux 



Figure 4. 138: Variation of Nusselt number with angle for Re=100,n=0.8 and Pr=l 

for constant heat flux 
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Figure 4.139:Variation of Nusselt number with angle for Re=10,n=0.6 and Pr=l 

for constant heat flux 



Figure 4. 140: Variation of Nusselt number with angle for Re=10,n=0.6 and Pr=10 

for constant heat flux 
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Figure 4.141:Variation of Nusselt number with angle for Re=10,n=0.6 and Pr=100 

for constant heat flux 
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Figure 4. 142; Variation of Nusselt number with angle for Re=10,n=0.6 and Pr=500 

for constant heat flux 







Figure 4. 143: Variation of Nusselt number with angle for Re=100,n=0.6 and Pr=l 

for constant heat flux 



Figure 4. 144; Variation of Nusselt number with angle for Re=100,n=0.6 and Pr=10 

for constant heat flux 
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Figure 4. 145; Variation ofNusselt number with angle for Re=l00,n=0.6 and Pr=100 

for constant heat flux 






Figure 4. 146: Variation of Nusselt number with angle for Re=100,n=0.6 and Pr=500 

for constant heat flux 
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Figure 4. 147: Variation of Nusselt number with angle for Re=l,n=0.5 and Pr=l 

for constant heat flux 
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Figure 4.148:Variation of Nusselt number with angle for Re=l,n=0.5 and Pr=10 

for constant heat flux 
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Figure 4. 153: Variation of Nusselt number with angle for Re=100,n=0.5 and Pr=100 

for constant heat flux 



Figure 4. 154: Variation of Nusselt number with angle for Re=100,n=0.5 and Pr=500 

for constant heat flux 


129 







Figure 4.158. 


Isotherms for Re=1 .0,n=1 .0.and porosity 0.4 Pr=1 .0(const wall temperature) 



Figure 4.159. 

Isotherms for Re=1 .0,n=1 .0.and porosity 0.4 Pr=10.0(const wall temperature) 
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Figure 4.160 



Figure 4.161 

Isotherms for Re=1 .0,n=1 .0.and porosity 0.4 Pr=1 00.0(const wall ten^rature) 
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isotherms for Re=10.0,n=1 .0, 
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Figure 4.168 


Isotherms for Re=100.0.n-1 .0.ar>d porosity 0.4 Pr-t .Ofconst wall temperature) 



Figure 4.169 


Isotherms for R0a1OO.O,n«1 .0.and porosity 0.4 Pral0.0(con8t wall temperatyr®) 
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Figure 4.180 

isotherms for Re=500.0,n=1 .0.and porosity 0.4 Pr-50.0(const wall temperature) 



Figure 4.181 


Isotherms for Rea500.0,n»1 .0.and porosity 0.4 Pr*1(}0.0(const wall temperatufi) 



Figure 4.182 

Isotherms for Re=1 .0,n=0.8,and porosity 0.4 Pn=1 .0(const wall temperature) 



Figure 4.183 


Isotherms for Re=1.0,n=0.8.and porosity 0.4 Pr=10.0(constwaii temperature) 
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Figure 4.194 


Isotherms for Re=1 00.0, n=0.8, and porosity 0.4 Pr=50.0(const wall temperature) 



Figure 4.195 


Isotherms for Re=1 00.0, nsO.8, and porosity 0.4 Pr=100.0(const wall temperature) 
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Figure 4.217 

Isotherms for Re-200.0, n-0.6.and porosity 0.4 Pr«10.0(cont wa« temperature) 













Figure 4.228 

Isotherms for Re«1 .0.n«0.5.and porosity 0.4 Pr-50.0{(»n8t wall temperature) 



Figure 4.229 


Isotherms for Re«1.0»r>«0.S,ahd porosity 0.4 Pr«® lOO.OCotmsi wall ttnrpifityr®) 
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Figure 4.236 

Isotherms for Res 1 00.0, n»0.5.and porosity 0.4 Pr»1 .0(cx>nst waff tefTperature| 





Figure 4.241 


Isotherms for Re=1 .0,n=0.8,and porosity 0.5 Pr=1.0(const wall temperature) 



Figure 4.242 

Isotherms for Rest. 0, 0*0.8, and porosity 0.5 Pra10.0(const wall temperature) 
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Figure 4.243 


Isotherms for Re=1 .0,n»0.8,and porosity 0.5 Pr»S0.0(conBt wall temperature) 



Figure 4.244 

Isotherms for Re-1.0,n»0.8,and porosity 0.5 Pr»100.0(oonst wall tsn^erature) 
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Figure 4.245 


Isotherms for Re=1 .0,n=0.8,ancl porosity 0.5 Pr=100.0(const wall temperature) 



Figure 4.246 


isotherms for Real .0, 0=0.8, and porosity 0.5 Pr=500.0(const wall temperature) 
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•e 4.262 
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Figure 4.269 


Isotherms for Re=1 .0,n^0.6,ancl porosity 0.5 Pr=100.0(const wall temperature) 
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Figure 4.273 


Isotherms for Re=1 0.0, 0=0.6, and porosity 0.5 Pr=50.0(const wall temperature) 



Figure 4.274 

Isotherms for Re=10.0,n=0.6,and porosity 0.5 Pr=100.0(const wall temperature) 
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Figure 4.275 


Isotherms for Re-10.0, n«0.6.anci porosity 0.5 Pr-500.0{const wall temperatyre) 



Figure 4.276 

Isotlienns for R®«*100.0,ii«0.63ocf porosity 05 Pf*1 .0fcon.st wilt 
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Figure 4.288 


Isotherms for R6=1.0,n=0.5,anc( porosity 0.5 Pr=:1 .0(const wall temperature) 



Figure 4.289 

Isotherms for Re=1 .0,n=*0.5,and porosity 0.5 Pr=10.0(const wall temperature) 
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Figure 4.292 

Isotherms for Re=1.0, 0=0.5, and porosity 0.5 Pr=500.0(const wall temperature) 








Figure 4.298 


Isotherms for Re»100.n=0.5,ancl porosity 0,6 Pt«i Ofwnt wall inmp) 



Figure 4.299 


Isotherms for Re«100,n«0.5,ancl porosity 05 Pr« lO Ofront well temp) 
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Figure 4.304 


Isotherms for R6=200.0,ns:0.6,and porosity 0.5 Pr=10.0(const wall temperature) 





Figure 4.305 

Isotherms for Re«200.0,n»0.5,ancl porosity 0.5 Pr*50.0(const wall temperature) 
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Figure 4.306 


Isotherms for Re«200.0.n.0.5.and porosity O.S Pr«1 00 Ofconst wall K'mpnf.it.jrf, 


Figure 4.307 


I'SOthemna for R@«200.0,n«0.S,af^d :pof»ity 0 J Pf»-SC)0.0fa»st w®ll liirptfifiiril 


\ 


j 
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Figure 4.312 


Isotherms for Re=1 .0,n=1 .0.and porosity 0.6 Pr=1.0(const wall temperature) 



Figure 4.313 

Isotherms for Re=:1 ,0,n=1 .0.and porosity 0.6 Pr=10.0(const wall temperature) 



Figure 4.314 

isotherms for Ra-1.0.n.1.0.andpomrty0.6Pr=^.0(cons^ 
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Figure 4.316 


Isotherms for Re»1 .0,n»1.0,ancl porosity 0.6 Prs500.0(cx)nst wall temperature) 



Figure 4.317 


Isotherms for Re=10.0,n=:1.0,and porosity 0.6 Pr=:1 .0(const wall temperature) 
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Figure 4.320 


Isotherms for R©»10.0,n»1.0,and porosity 0.6 Pr^100.0(const wall temperature) 



Figure 4.321 

Isotherms for Re=10.0,n=1 .0.and porosity 0.6 Pr=500.0(const wall temperature) 
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Figure 4.322 

Isotherms for Re-IOO.O.n-l .O.and porosity 0.6 Pr-1 .0{con8t wall temperature) 
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Figure 4.328 


Isotherms for Re=200.0,na1.0,and porosity 0.6 Pr=!l0.0(const wail temperature) 



Figure 4.329 

Isotherms for Re=:200.0,n=1 .0,ancl porosity 0.6 Pr=50.0(const wall temperature) 



217 



Figure 4.330 


Isothems (or Re.200.0.n-1.0.and poroaity 0.6 Pr-100 0(c»n« wafl 



Figure 4.331 


tsothermsf for Rt«200.0,n*«1 .O^aod O.S Pr*5W.0|oonil will 


Figure 4.332 


Isotherms for Re«=500.0,n=1.0,and porosity 0.6 Pr=1.0(const wall temperature) 
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Figure 4.334 


isotherms for Re-500.0, n-1 .0.snd porosity 0.6 


pr!.5O.0(oonst wall tompotalurs) 



Figure 4.335 


isotherms for Re-500, 0.n-I.O, and pofwiity 0.6 Pr-100 Ofoonst wa« tempiwatureS 
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Figure 4.336 


Isotherms for R6=1 .0,n=0.8,and porosity 0,6,Pr=1 .0(const wall temperature) 



Figure 4.337 


Isotherms for Real .0,n=0.8,and porosity 0.6.Praio.O(const wall temperature) 




Figure 4.338 


isotherms for Re-1 .0.n.0.8,and porosity 0.6.Pr.50.0(cor,st wall temperature) 



Figure 4.339 


isotherms for Re-IO.O.n-O.S.and porosity 0.6 Pr.100.0(«.hSt wall temperature) 
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Figure 4.344 

Isotherms for Re=»10.0,ns=0.8,and porosity 0.6 Pr*s100.0(c5onst wall temperature) 
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Figure 4.354 

isotherms for Re-200.0.rt=0.8,arKi porosity 0.6 Pr.lOO.Ofoonst w 



Figure 4.355 


Isotherms for Re-200.0, n-0.8.and porosity 0.6 Pr-SOO.Ofoonst 



Figure 4.356 

Isotherms for Re=t500.0,n=s0.8.afKJ porosity 0.6 Prx:1.0(const wall temperature) 



Figure 4.360 


Isotherms for Ress1 .0,n=0.6,and porosity 0.6 Pr=1 .0(cx)nst wall temperature) 



Figure 4.361 


Isotherms for Real .0,n=0.6,and porosity 0.6 Praio.0(const wall temperature) 



Figure 4.362 


Isotherms (or Re-1 .0.n-O.e.arxJ porosity 0.6 Pr»50.0(ooosl wall temperature) 



Figure 4.363 


Isotherms for Re*1 .0,tt«0.6,BrKj porosity 0.6 Pr»1 00.0(oon®t wall temperature) 
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Figure 4.364 


Isotherms for Re=1 .0,n=0.6,and porosity 0.6 Pr=500.0(const wall temperature) 



Figure 4.365 


Isotherms for Re=10.0,n=0.6,and porosity 0.6 Pr=1 .0(const wall temperature) 
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Figure 4.366 


Isotherms lor Re«10,0,n»0.6.and porosity 0 6 Pr«t0 Ofconst wall tetrrperaturo) 
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Figure 4.368 

Isotherms for Re=10.0,n=0.6,and porosity 0.6 Pr=1 00.0(const wall temperature) 
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Figure 4.372 


Isotherms for Re=1 00.0, n=0.6, and porosity 0.6 Pr=50.0(const wall temperature) 





Figure 4.373 

Isotherms for Re«T00.0,n=0.6,and porosity 0.6 Pr«100.0(const wall temperature) 
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Figure 4.374 


Isotherms for Re=100.0,n^.6.and porosity 0.6 Pr.SOO Ofconst wall tomppratyrei 



Figure 4.375 

Isotherms for Rea200.0,na0.6,afKl porostty 0.6 Pr*1 Ofc^ormf wall mmmt Blum) 


Figure 4.376 


Isotherms for Re=200.0,n=0.6,and porosity 0.6 Pr=10.0(const wall temperature) 



Figure 4.377 

Isotherms for Re=200.0,n=0.6,and porosity 0.6 Pr=50.0(const wall temperature) 
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Figure 4.378 


Isotherms for Re*200.0,n=0.6,and porosity 0.6 Pr«100 Ofconst wall temperature} 



Figure 4.379 


Isotherms for Rea2O0.0*n*0vS,aoc( .porosity 0.6 Pfs5CX).0|a)«st wait 


Figure 4.380 

Isotherms for Re=:500.0,n=0.6,and porosity 0.6 Pr=1.0(const wall temperature) 



Figure 4.384 


Isotherms for Re=1 .0, 0=0.5, and porosity 0.6 Pr=1.0(const wall temperature) 



Figure 4.385 

Isotherms for Real .0,n=0.5,and porosity 0.6 Pr=10.0(const wall temperature) 




246 


Figure 4.388 


Isotherms for Re=1 .0, 0=0.5, and porosity 0.6 Pr=500.0(const wall temperature) 



Figure 4.389 


Isotherms for Re=10 0,n=0.5,and porosity 0.6 Pr=1.0(const wall temperature) 
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Figure 4.392 


Isotherms for Re=a10.0,n=0.5,and porosity 0.6 Pr=:100.0(const wall temperature) 
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Figure 4.396 


Isotherms for Re*s1 00,0, na=0.5, and porosity 0.6 Pr=50.0{const wall temperatura) 





Figure 4.397 

Isotherms for Re*! 00.0, naO.5, and porosity 0.6 Pr=100.0(const wall temperature) 
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Figure 4.412 

Isotherms tor Re-100 0,n»l O.and porosity 0 ■» Pr»it» host fi,j. 



Figure 4.413 


Isothtrms for :Rt»2(KJ.0sn*1..0,itin 0.4 P'r»i0.Clfi«fml fttil «ii« I 
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Figure 4.428 


fm R®«iO O.nwO 8.*mi pomiiuy 0.4 Pf»lO Of'mmt »%»»f ftw»| 



Figure 4.429 


Isothiorma for Pi;«««l0,0.rt«n>.0.#,arn!J poroaily 0.4. .F»f«.,SO.O|«:»ttl I1M*I 
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Figni'P 4.448 

Isotherms for Re»iO 0.n»0.6,and porosity 0 4 Fr-sr, n.roost heal 



Figure 4.449 


IsoUiirms for Rt^10.0.n4l.6,anc3 pofowlf 0.4 fmm 




Isotherms for Re«100,0.ns0 6,andporosilv0 4 Pr«!'!(ir’'«iir<ii <( .» 



Figure 4.453 

Isotherms lot Re* tot) 0,n»0 B.at\rt porostty 0 4 Pr«50 Ofronst h»ii flu* I 
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Figure 1.460 

Isotherms (or Re*200 0.t>:=0 6,atitl poros.fV' 0 4 P- ■■ 



.‘•,r e*.* Fr.l* * ..!l 
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Figure 4.465 


Isotherms for Re=1 ,0,n=0.5,and porosity 0.4 Pr=500.0(const heat flux) 



Figure 4.466 


Isotherms for Rej!l0.0,n=0.5,and porosity 0.4 Pr=100.0(corist hsat flux) 
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Fii-iu. f !' 





Figure ‘1.408 

Isotherms (or Re* 10 0.n«0 S.and potowty n 4 Pt-.'iixj tiiijunsi heal itu* ) 
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Figure 4.469 

Isotherms for Re=1 00.0, n-D.S, and porosity 0.4 Pr=1 .0(const heat flux) 



Figure 4.470 


Isotherms for Res 1 00.0, nsO.S, and porosity 0,4 Pr=10.0(const heat flux) 












isotherms tor Re-IOO.O.n^ S.end porosrty 0 5 Pr- tOO Olarnt ho.! «u.l 



Figure 4.489 

Isothertnt lor R*»100 0,ft»0 B.end powity 0 5 Pr»100 Ojoooirt heel «o»l 
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Figure 4.512 

Isothtrms fw l»rosiff 0.S litil ii«| 



Figure 4 J13 


iiothtrms for Fi(i*2CK),fi»0.®, tr^ peioiitf §.| ef'«10 CMcoti w»l t«if iiwl 




311 














4.530 


iv.OwTm-, f.,ft.!;r,,j)nflf)MWily05Pf.1!»,0(conslheatnux) 
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Figure 4,533 


-n. ir, r,*. - w, n-0 i 00,^ 0 S ocoort m tmt flux) 



Figure 4.534 


iir«h*r«^i* lao jtwO 0 5 Pr»100 0*(W« w»8 
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Figure 4.537 


Isotherms for Re=200.0,n=0.5,and porosity 0.5 Pr=1 0.0(oonst heat flux) 



Figure 4.538 

Isotherms for Re«200.0,h«0.5,and porosity 0.5 Pr=50.0(const heat flux) 
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Figure 4.539 

Isotherms for Re»200.0,n=0.5.and porosity 0.6 Pr«l 00.0(oon8t heat flux) 



Figure 4.540 


Isothemw for Re«200.0,n«0.5.«f)d porosity 0.6 Pr«600.0(oonst hMt flux) 
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Figure 4.557 


Isotherms for Re=100.0,na1 .0,ancl porosity 0.6 Pr=50.0(const heat flux) 



Figure 4.558 

Isotherms for Re«100.0,n=1 .0.and porosity 0.6 Pna100,0(const heat flux) 



333 




Figure 4.558 


for ^|>ort»fty 0 J Pr«1 .OCcwnM Mm MkI 


Figure 4.559 


Isotherms for Re*=200.0,n*1.0,ancl porosity 0.6 Pri=10.0(con8t heat flux) 



Figure 4.560 

IsotheiTTW for Re«200.0,n«1 .0,and porosity 0.6 Prte50.0(c»n8t heat flux) 
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Figure 4.620 


lor (^100 e.rMJ pr»rM«v 0 « ^»»Vt twin »j,» ) 



Figure 4.621 


iK^wftw tor II|N100,0/I-O.8jmtf pMMy 0 B Pf.iOO OfmM Nkm But} 




Figure 4.624 


^ © f ^*1©,§|»iil *i»;| 



Figure 4.625 

!»o«»mw Jw I^Ni5SJO.O,iM).8.*nd ptxogty OJ 0(«rai (wai Ami) 
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Figure 4.628 

Isothemw for Re- 500 . 0 , M).5,and porosily 0.6 Pf»1 OfoonsJ heat fiu*) 



Figure 4.629 


leothenw tor R»-500.0,n-0.5jtftd porcwHy 0.6 Pr-l0.0{<»n« f»**l Hu*) 



Figure 4.630 

*r» «|»*900 S.ml poRMty 0.8 I^.Ofoonw h»al to) 



Figure 4.631 


iwtwfw *m l%»a00.O.iW».#,i«l pwo^ 0,8 P»»i100.0(axi8f heat flux) 
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Chapter 5 


CONCLUSIONS AND SUGGESTIONS 
FOR FUTURE WORK 


5.1 Conclusions 


In this work, the incompressible steady flow and heat transfer of power la\N 
fluid across a bank of circular cylinders has been studied numerically. The inter- 
cylinder interactions have been simulated using a simple cell model which does noi 
exchange mass or energy with the neighboring cells. The resulting equations oJ 
continuity, momentum and energy have been solved using the fini te differenc€ 
method to obtain the velocity and temperature fields. These in turn have beer 
further processed numerically to obtain the values of the local and average Nusseli 
number over a wide range of physical and kinematic parameters as : 
0.5s n S 1. 0,0.40s e S 0.60,1 S Re S 500 and iSPrSSOO for both boundar> 
conditions, i.e., the constant wall temperature and constant heat flux at the surface oJ 
the cylinder. The solution procedure and computer code have been validated usin^ 
the results for limiting conditioiK available in literature for both Newtonian anc 
non-Newtonian fluids. 

The observed dependencies of the Nusselt number on porosity, Reynolds 
number, Prandtl number , power law index and boundary conditions have beer 
explained physically. As expected, the value of the Nusselt number is always higher ir 
the case of constant heat flux condition than that for the isothermal case, though the 
extent of this enhancement is strongly dependent on the value of the porosity, 
Reynolds number , power law index and Prandtl number. At low Peclet numbers, the 
temperature field shows complete fore and aft symmetry which is progressively losi 
with the rising value of the Reynolds number and / or Prandtl number which indicates 
the increasing role of convective contribution to heat transfer. 

The Nusselt number is influenced by the power law index, i.e. shearing 
action of fluid. As the value of tiie power law index decreases, the rate of heat trmisfei 
increases, because shear thinning lowers the resistance to heat flow and this facilitates 
heat transfer. Also the results, at low values of Reynolds number, are not muct 
influenced by power law index. However, as the value of the Reynolds numbei 
progressively rises, shearing action influences the heat transfer rate. 
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1 he value 0i the Nusscit number is also strongly influenced by the value of 
porosity. The dependence is consistent with flie intuitive feeling that owing to steeper 
velocity and tc?nperaturc gradients in concentrated ( low porosity) systems, the value 
of the Kusselt number shows an inverse dependence on the porosity. The present 
theoretical results have been validated using the scant experimental data available' in 
the literature. 

5.2 Suggestions for future work 

Considerable scope exists to extend tJus work. Some of the possibilities are 

mentioned here. 

1 . Although the power law fluid model gives the simplest representation of 
shear-thinning non-Newtonian behavior, its inability to predict a constant 
value of viscosity at the outer surface, cast doubt about its appropriateness 
to describe flows with stagnation points. Thus one must use a more 
realistic model such as the Ellis or the Carreau equation. 

2. Similar computations can be performed for more realistic non-isothermal 
processes of non-Newtonian fluids heatbansfer from the. cylinders. 



BIBLIOGRAPHY 


Adams, Don., and Bell, K.J., Fluid Friction and Heat Transfer For Flow of Sodium 
Carboxymethylcellulose Solutions Across Banks of Tubes, Chem. Eng. Prog. Sym. 
Ser. ,No. 82, Vol 64, 1968, pp.133-145 

Ahmad, R.A., Heat Transfer Engineering, Vol 17, 1996, pp.31-81 


Bailer, F., Tochon, P., Grilliot, J.M. and Mercier, P., Rev. Gen. Therm., Vol 36, 
1997, pp-744-754 

Dmmmmond, J. E., and Tahir,M.I., Int. Jl. Multiphase Flow, 10,515, (1984) 

Happel J., Viscous flow relative to array of cylinders. AIChE J. 5,p. 174, (1959). 

Hsu, C.- J., International jounal of Heat and Mass Transfer, Vol 7, 1964, pp. 431- 
446 

Kreith, F. (ed), The CRC Handbook of Thermal Engineering, CRC Press, Boca 
Raton, FL(U.S.A.) pp.3.26-3.46 

Lange, C.F., Durst, F. and Breuer, M., International Journal of Heat and Mass 
Transfer, Vol 41, 1998, pp.3409-3430 

Mandhani, V.K.,Chabbra, R.P and Eswaran, V., Forced Convection Heat Transfer 
in Tube Banks in Cross-Flow, Chem. Engng. Set, in press (2002). 

Martin, A.R., saltiel, C. and Shyy, W., International Journal of Heat and Mass 
Transfer, Vol 41, 1998, pp.2383-2397 

Met 2 ner, A. B., and J.C. Reed, AIChE J., Vol 1, 1955, pp.434-440 

Metzner, A. B.,R. D. Vaughn, and G. L. Houghton., AIChE J., Vol 3, 1957, pp.92- 
100 

Nishimura, T., Itoh, H. and Miyashita, H., International Journal of Heat and Mass 
Transfer, Vol 36, 1993, pp.553-563 

Sangani, A. and Acrivos, A., Saltiel, C. and Shyy, W. , International Journal of 
Heat and Mass Transfer, Vol 41, 1998, pp.3409-3430. 


374 



Shah, M. N., E. E. Petersen, and Andreas Acrivos., AIChE J., Vol 8, 1962, pp.542- 
549 

Skartsis, L., Khomammi, B., and Kardos, J.L.//. Polymer. Eng. Set, 23, 32, 
(1992) 

Tripathi, A., and Chhabra, R.P., Transverse laminar flow of non-Newtonian fluids 
over a bank of long cylinders, Chem. Engg. Commm. Vol. 147, pp 197-212, 
(1996) 

Tripathi, A., and Chhabra, R.P., Slow power law fluid flow relative to an array of 
infinite cylinders, Ind, Eng. Chem. Res. Vol 3 1 , pp 2754-2759, (1992) 

Zukaukas, A., Convective Heat Transfer in Cross-Flow in Handbook of Single- 

Phase Convective Heat Transfer (edited by S.Kakac.R.K.Shah, W.Aung) Wiley- 

Interscience , NewYork(1987). 


375 



APPENDIX - A 


Adams and Bell(1968) work-relating the j ^ factor with present parameters 
for a special cases of porosity 0.714,n=0.83 and Pr=740.56 (0.5% CMC 
solution) and porosity 0.75, n=0.748 and Pi=9921.36(l%CMC solution). 

Fluid fiiction and heat transfer data are presented for flow of dilute 
sodium carbox3miethylcellulose (CMC) solutions in water across the ideal 
tube banks containing an equilateral triangular tube layout and in-line 
square arrangement. The tube banks contained 0.375in. O.D. tubes in a 1.25 
pitch ratio. 

The Reynolds number used, for the calculation part, is 


DT^’''p 


Re 


A-B 


r 


where F = 


U 


(1-Vi^) 




^ 4n y 


(A.1) 

(A.2) 

(A.3) 


relating with the present Reynolds number, Re = 


m 


Re^.j =ReF(n,£) 
where, F(n,s)= — 


1 


8 


n-l 


3w + 1 Y 
4n J 


(A.4) 

(A.5) 


(l-VT^)^"" 


The Prandtl number used by Adams and Bell is 


^^A-B 


K [D 


(A.6) 


C m(uy ^ 

relating with the present Prandtl number, Pr = • 


JC 




Pr . .=PrF,(n,£) 


Where, F, («,£) = 8 


n-l 



y, f \ 

1 1 

4n ) 

li-vr=Tj 


n-l 


(A.7) 

(A.8) 
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The heat transfer results are presented in plots of factor vs. Reed- 
Metzner Reynolds number as. 


Jh 


C,>GA^ L 


£^(L 
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. 11-1 


2/3 


(A.9) 


where _ ^ 3/? + 1 _ non - NewtonianfluidNusseltnumber 

1 4/2 ' 'KT^'y £ ft -i ^ , 2 ^]\T^ . ^ ^ J. 7_ ^ 


NewtonianfluidNusseltnumber 


10 ) 


Relating the factor factor with Re ,Nu and Pr; 


Using the difmation of Nusselt number, Nu= — 

k 

Re and Pr , the j„ factor can be written as 
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where, F 2 (n,e) = S 
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(A. 12) 

(A.13) 

(A. 14) 

(A.15) 

(A.16) 
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A(a). Model 1 of the Adams and Bell (1968): 

For the special case of porosity of 0.75(n=0.748) with the tube banks of an 
equilateral triangular arrangement and for 1% CMC solution(Pr=921.36) , 
the above relations becomes 


Re^_a=3.78616Re 

Pr,.^=0.5282396Pr 

(A. 17) 

(A. 18) 

=0.31804182 

RePr*'^ 

The experimental j„ factor can be related by 

(A.19) 

j„ = 0.765 Re"°“’- 0.0086 

(A.20) 

A(b). Model 2 of the Adams and Bell (1968): 


For the special case of porosity of 0.712(n=0.83) with the tube banks of 
in-line square arrangement and for 0.5% CMC solution(Pr=740.56) , the 
above relations becomes 
the relations becomes 

Re =3.3448 Re 

Pr^.^ =0.64265 Pr 

(A.21) 

(A.22) 

/„ = 0.22253 

RePr‘^" 

The experimental factor can be related by 

(A.23) 

7^ =0.599 Re^-®*’+ 0.003 

(A.24) 
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